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1. Introduction 

In random matrix theory, the determinant is naturally an important functional. The study 
of determinants of random matrices has a long history. The earlier papers focused on the 
determinant detv4 n of a non-Hermitian iid matrix A n , where the entries of the matrix were 
independent random variables with mean and variance 1. Szekeres and Turan [23] studied an 
extremal problem. Later, in a series of papers moments of the determinants were computed, see 
[20J and [4] and references therein. In [24] , Tao and Vu proved for Bernoulli random matrices, 
that with probability tending to one as n tends to infinity 

\/n7exp(— cyn logn) < |detA n | < \fn\ u(n) (1.1) 

for any function u(n) tending to infinity with n. This shows that almost surely, log | det A n \ is 
(~ + o(l))n logn. In Goodman considered the random Gaussian case, where the entries 
of A n are iid standard real Gaussian variables. Here the square of the determinant can be 
expressed as a product of independent chi-square variables and it was proved that 

log (| det A n \) - 7, logn! + f logn 



\ logn 



->iV(0,l) R , (1.2) 



where N(0, 1)r denotes the real standard Gaussian (convergence in distribution). A similar 
analysis also works for complex Gaussian matrices, in which the entries remain jointly indepen- 
dent but now have the distribution of the complex Gaussian N(0, l)c- In this case a slightly 
different law holds true: 

logfl det A n \) - \ \ogn\ + \ \ogn 



i logn 



->AT(0,1) R . (1.3) 



Girko [H] stated that (II. 2p holds for real iid matrices under the assumption that the fourth 
moment of the atom variables is 3. In [10] he claimed the same result under the assumption 
that the atom variables have bounded (4 + 5)-th moment. Recently, Nguyen and Vu [19] gave 
a proof for (11.21) under an exponential decay hypothesis on the entries. They also present an 
estimate for the rate of convergence, which is that the Kolmogorov distance of the distribution 
of the left hand side of (II. 2p and the standard real Gaussian can be bounded by log~3 +0 W n . 
In our paper we will be able to improve the bound to log" 3 n in the Gaussian case. 

In the non-Hermitian iid model A n it is a crucial fact that the rows of the matrix are jointly 
independent. This independence no longer holds true for Hermitian random matrices, which 
makes the analysis of determinants of Hermitian random matrices more challenging. The 
analogue of (II. ip for Hermitian random matrices was first proved in [251 Theorem 31] as a 
consequence of the famous Four Moment Theorem. Even in the Gaussian case, it is not simple 
to prove an analogue of the Central Limit Theorem (CLT) (II. 3p . The observations in [11] do 
not apply due to the dependence between the rows. In [18] and in [T4], the authors computed 
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the moment generating function of the log-determinant for the Gaussian unitary and Gaussian 
orthogonal ensembles, respectively, and discussed the central limit theorem via the method of 
cumulants (see [TJJ equation (40) and Appendix D]): consider a Hermitian n x n matrix X n in 
which the atom distribution Qj are given by the complex Gaussian N(0, l)c for i < j and the 
real Gaussian iV(0, 1)r for i = j (which is called the Gaussian Unitary Ensemble (GUE)). The 
calculations in [H] should imply a Central Limit Theorem (see Remark 12.41 in our paper): 



Recently, Tao and Vu [26J presented a different approach to prove this result approximating 
the log-determinant as a sum of weakly dependent terms, based on analyzing a tridiagonal 
form of the GUE due to Trotter [27]. They have to apply stochastic calculus and a martingale 
central limit theorem to get their result. This method is quite different and also quite involved. 
More important for us, the techniques due to Tao and Vu seem not to be applicable to get 
finer asymptotics like Cramer-type moderate deviations, Berry-Esseen bounds and moderate 
deviations principles. The reason for this is the quality of the approximation by a sum of 
weakly dependent terms they have chosen is not sharp enough. Let us emphasize that Tao and 
Vu proved the CLT (11. 4ft for certain Wigner matrices, generating a Four Moment Theorem for 
determinants. 

The aim of our paper is to use a closed formula for the moments of the determinant of a GUE 
matrix, giving at the same time a closed formula for the cumulant generating function of the 
log-determinant. We will be able to present good bounds for all cumulants. As a consequence 
we will obtain Cramer-type moderate deviations, Berry-Esseen bounds and moderate deviation 
principle (for definitions see Section 2) for the log-determinant of the GUE, improving results 
in [H] and [26]. Moreover we will obtain similar results for the GOE ensemble. Good estimates 
on the cumulants imply such results. To do so we apply a celebrated lemma of the theory 
of large deviations probabilities due to Rudzkis, Saulis and Statulevicius [21], [22] as well as 
results on moderate deviation principles via cumulants due to the authors [6]. Applying the 
recent Four Moment theorem for determinants due to Tao and Vu [26], we are able to prove 
the moderate deviation principle and Berry-Esseen bounds for the log-determinant for Wigner 
matrices matching four moments with either the GUE or GOE ensemble. Moreover we will be 
able to prove moderate deviations results and will improve the Berry-Esseen type bounds in [19] 
in the cases of non-symmetric and non-Hermitian Gaussian random matrices, called Ginibre 
ensembles. 

Remark that the first universal result of a moderate deviations principle was proved in [8] 
and [7] for the number of eigenvalues of a Wigner matrix, based on fine asymptotics of the 



log( | det X n \ ) - | log n! + \ log n 




(1.4) 
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variance of the eigenvalue counting function of GUE matrices, on the Four Moment theorem 
and on localization results. 

2. Gaussian ensembles and Wigner matrices 

Among the ensembles of n x n random matrices X n , Gaussian orthogonal and unitary ensem- 
bles have been studied extensively and are still being investigated. Their probability densities 
are proportional to exp(— tr(X^)), where tr denotes the trace. Matrices are real symmetric for 
the Gaussian orthogonal ensemble (GOE) and Hermitian for the Gaussian unitary ensemble 
(GUE). The joint distributions of eigenvalues for the Gaussian ensembles are ([U Theorem 
2.5.2], [13 Chapter 3]) 



n ' p v i=l 7 l<j<k<n 



P n , /3 (A 1 ,...,A n ):=^exp(-^>^A l 2 ) [\ \Xj — \ k \P, (2.1) 

l<j< 



where /9 = 1, 2 for the orthogonal and unitary ensembles, respectively, and Z n ^ is the normaliz- 
ing constant, sometimes called the Mehta integral (see [H Theorem 2.5.2, formula (2.5.4), and 
Corollay 2.5.9, Selberg's integral formula]). 

Let us denote by X@ the random matrices of the two Gaussian ensembles. We are interested 
in the moments of | det for these ensembles, that is 

„ n 

M^(s):=(|detXf|V= / *W>(Ai, ■ ■ ■ , A n ) TT |A,| S d^. 

All information about the distribution of log|detX^| can be obtained from the generating 
function M n $(s). The moments of log|detX^| may be obtained from the coefficients in the 
Taylor expansion of M n>i g evaluated at s = 0, 



((logldetXfl)^ 



j>0 



the corresponding cumulants Tj(n, (3) := (— i) J 'Jp- logE[V' log ' detX ™'] | are related to the Tay- 
lor coefficients of log M n ^ via 



log M ntP (s) =J2 



r^) 



v- 



In the literature the Mellin transform of the probability density of |detX^| was calculated 
for the Gaussian ensembles, giving an explicit formula for M n ^(s). To be more precise, if 
9n,&{') denotes the probability density of the determinant of a GOE or a GUE matrix and 
dnsiy) := \{9n,p{y) + 9n,p(—y)) be the even part, the Mellin transform of g^a is defined by 



oo 

8-1 



MnAs) ■= / y s 9nAy)dy- 
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For the GOE and GUE ensembles we obtain 

^ rco PCO 

M n ,p(s) = - •••/ P n: fs(X 1 ,...,\ ri )\X 1 ---X n \ s - l d\ 1 ---dX n 

* J — CO J —CO 

and an obvious consequence is the relation 

M n ^s) = 2M n ^s + l). (2.2) 

It is quite involved to calculate the Mellin transform even for the Gaussian ensembles. The case 
— 1 was calculated in [HI formulas (31), (19) and (26)] (see also [UJ Chapter 26.5]). Here 
the Mellin transform is a Pfaffian of an anti-symmetric matrix applying the method of (skew) 
orthogonal polynomials. With (12. 2p . for n = 2p + 1 one obtains 

M w {s) = 4-/2 Yl [ ^l + hl f (2.3) 

m=l 1 U + °m) 

with b x m := | + z- If n — 2p one obtains 

^ p fs + l s n + l + s l^ r((s + l)/2)r((n+l)/2) Tjs + m + p 

M ^ s) ~ 2 F y—-*—2—^2) r(i)r((n + i + s) /2) 11 i> + i) ' 

(2.4) 

where F is the (GauB) hypergeometric function 

m=0 w 

with (a;)*" 1 ) := x(x+l)(x+2) ■ ■ ■ (x+m — 1) denoting the Pochhammer symbol. F is convergent 
for arbitrary a, b, c and for real — 1 < z < 1. In [3], an alternative derivation for (12. 3 p and (12. 4p 
is presented using terminating hypergeometric series. The case — 2 was calculated in [18, 
Section 2]. Here a knowledge of determinants and orthogonal polynomials is needed. One 
obtains 

M n , 2 (s) = 2-/ 2 J] [ ^l + h2) m) (2.6) 

m=l ^2 "I" "mi 

with 6^ = LyJ • As a consequence of (12. 6 j) we obtain the following results for the cumulants 
1^,2) of log | det^|: 

Lemma 2.1 (Bounds for the cumulants of log | det GUE). For the Gaussian unitary en- 
semble = 2 we obtain 

r\(n,2) = --(l + l og 2) + -log(2Ln/2j) + const +0(1 /n) 

and 

a\ := r 2 (n, 2) = ~ log(2Ln/2j) + ~( 7 + log 2 + 1) + 0(l/n), 
where 7 denotes the Euler-Mascheroni constant. Moreover for any j > 3 we have 

|rj(n,2)| < const j\. (2.7) 
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Proof. Let us remark that some of our calculations can be found in [Mj. We work out all the 
details to get good bounds on the cumulants, which is not the aim in [13]. With ip(x) := 
^logr(x) we denote the dig amma function. From ( 12. 6 p we obtain 



r!(n,2) = ^logM ni2 (s) 



1 n 

log2 + -£^(1/2 + 6?). (2.1 



n 2 ° 2 
s=0 j = i 



For any n = 2k + 1 we obtain § Y*=i ^(V 2 + & ?) = Ej=i ^(V 2 + j) + IV'd) and for n = 2A; 
we have 1^=1^(1/2 + ^) = £* =1 ^(1/2 +j) + §V>(V 2 ) - ^H 2 ^)- WithT(l + x) = xT(x) it 

follows that + = ip(x) + ^ and therefore recursively ip(l/2+j) = ip(l/2) + 2 ^X/f=i 2l^T^ ' 
see [HH Section 1.3, (1.3.9)]. Using 

j=l 2=1 2=1 2=1 v 2=l 2=1 7 

we obtain £* =1 ^(1/2 + j) = ^(1/2) - fc + (2* + 1) fe?=i f - E?=i i) • Applying 

" 1 11 

£ 7+ l 0gn+ +0( ), (2 .9) 

2=1 

it follows that (2k + 1) (j^ti 1 - Eti = ( 2A: + + 21og2) + (2k + 1)1 log fc + 0(1). 

With ^(1 /2) = -2 log 2 - 7 we have 

X>(l/2 + •?') + ^(V2) = -* + (* + ~) log A; + O(-). (2.10) 

i=i 

In the case n = 2k we have to consider in addition the term l^(l/2 + k) = |logA; + 0(\). 
Summarizing we obtain for every n: 

ri(n,2) = --(log2 + l) + - log(2fc) + const +0(1/ n). 

From ( 12. 6 j) and ( 12. 8ft we obtain for n = 2A; + 1 



d 7 ' 

r i (n,2) = — logM n , 2 ( S ; 



iv°- 1) (l/2) + E ^ (i_1) (l/2 + (2.11) 



s =o 2^ 2^ 



with the polygamma function ip^ k \x) := J^r logr(x). For n = 2k one has to subtract from the 
right hand side the term l^ -1 ^). We remind the representation of T(x)~ 1 due to Weier- 
strass (see for example |T6l Section 1.3, (1.3.17)]): = xe lx YYk =l (^- + f )e~k. Differentiating 
— logT(x) leads to 

^ )= - 7 -^ + £G-^)=- 7+ £(^tt-^)- 

fc=l V 7 n=0 V 7 
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Therefore one obtains 

oo 1 

^(x)^-!^^^— T . (2.12) 

n=0 ' 

It follows that 

k k oo 

^-^1/2 + i) = (-^(i-i^^EET^riy 

i=l i=l m=i ^ ' 

(k k k oo 1 \ 

2 ££<5^TT>7 + 2 £ £ 7>TTT)jh :r ' + T - 
i=l m=i v 7 8=1 m=fc+l v 1 / / 

With 2 Eti E™=* (2^TT)J = Eli (2m+ 1 i )J -i " ELi (2^W we obtain 

T x = (-I)' (j - 1)! 2 j - 1 V 7 ^— — - (j - 1)! 2^ - (-1V (j - 1)! 2^ V 

m=0 v ' m=l v 2 

Further we get 



(2m + IV' ' 

m=fc+l v 7 



Hence using (PZU for we obtain 



fc A; 

J>^)(l/2 + z) = (-l)i (j - 1)! 2*- 1 E f2m * lV -i ~ (h 

i=l m=0 ^ 2 

oo 

+ (-iy(j-i)\23-\2k + i) J2 ro ■ nj - ( 2 - 13 ) 

' 2m + IP 

m=fe+l v 2 

In particular for j = 2, we have 

E (V2 + = 2(1 log(*) + 1(7 + 2 log(2))) - 1^(1/2) + \(2k + 1)^ (k + |) 
i=i 

= log(A;)+ 7 + 21og(2)-^ (1) (l/2) + 1 + 0(1). (2.14) 
With fl2TLj) we obtain forn = 2A; + 1 that 

A: oo 1 

r j(B , 2) = (-!)*(,- i)i S^jjzr + Hyw- 1)" (2* + i) £ p^j- 

m=0 v 7 m=k+l y ' 

The first term is -2^3 (j -l)^- 2 )(l)+0(l/k). The second term is 2~ J '(2A; + l)^ " 1) (| + ^ + l)- 
For n = 2k we have to subtract 2~ 3 'ip^~ 1 \^ + fc). Finally we will apply some bounds for the 
polygamma functions ip^\ Therefore we will apply the following integral-representation (see 
for example [161 Section 1.4, (1.4.12)]): 

V>(a;) = log(a;) - / e~ tx (tf(t) + 1 j tft with /(t) :=(!-! + — i— - ] -, t > 0. (2.15) 



2 t e* - 1 J t 
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Differentiating we see that for j > 1: 



(2.16) 



Notice that < (tf(t) + ±) < 1 for every t > 0; hence we obtain for every x > and every 
3 > 1: 



|V> (i) (z)| < j\x- j + j\x~ j ~ l . 



(2.17) 



Let us consider the variance a\ = ^(n, 2). With (I27T7D we have |</> (1) (l/2 + k)\ < (§ + A;)" 1 + 
(| + A;)" 2 . Hence we have cr| = § V (1) ( 1 / 2 + i) + |^(l/2) + and with (EHJ) we 



obtain 



1 



1 



a 2 = - log fc + -( 7 + 2 log 2 + 1) + 0(l/k). 
For j > 3 the cumulants can be bounded by: With (12.171) we obtain 



|r,(n,2)| < 



+ 



2 _J '(2fc + l)^ (i_1) (l/2 + fc + 1) 



2~ j i/) ( i~ 1 \l/2 + k) 



0(l/k) 



< 6(j — 1)! + const 



2J- 1 A;^" 2 



+ 



2J- 1 jfeJ- 



i < const (j - 1) 



Therefore the cumulants satisfy the stated bounds. 



□ 



With some more technical effort we obtain similar results for the Gaussian orthogonal en- 
sembles: 

Lemma 2.2 (Bounds for the cumulants of log | detX*|, GOE). For the orthogonal Gaussian 
ensemble ((3 = 1) we obtain 



n n 
r x (n, 1) = -log(2|ri/2j) - - + const +0(1/ 



and 



of := T 2 (ft, 1) = log(2Lft/2j) + 2 + 1-27^ + ^ + 0(1 /n) 
where K denotes Catalan's constant K = J2m=o (^m+i) 2 > an< ^ f or an y 3 — ^ 

|r_j(ft, 1)1 < const j\. 



7T 
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Proof. For /3 = 1 and n — 2k + 1, formula (12. 3p for the Mellin transform implies 

r l( „,i) = ^i og M„»U = ^io g (4) + ix:*(i + i[i^j + i) 

i=l 

= nlog(2) + ^^( i + ^) + 2^(4 + 2) 

i=0 

= nIog (2) + l*(|) + t(i*(| + ^) + i*(| + |)). 

i=l v 7 

The last transformation is useful since we are now able to apply Legendre's duplication formula 
T(z)T(z + 1/2) = 2 1 - 2z ^tvT(2z) (see for example [16, Section 1.2]). This implies 

^(z) + l -^{z + i) = ^(2z) - log(2). (2.18) 
With z = 3/4 + i/2 - 1/2 we obtain 

fc 

r 1 (n,l)=nbg(2) + ^(f)+X)V'(V2 + i)- felo g(2)- (2-19) 

i=i 

The summand §"0(1) equals via the same identity t/>(§) — log(2) — \ip(j) = ^(f) — l°g(2) + 
f + §log(2) + §7 = f - |log(2) - ± 7 . As in the GUE case, we have Eti^(V 2 + *) = 
- k + (k + |) log(fc) + O(i), see fl2TTUj) . Now (gjg]) implies that 

ri(n, 1) = 2 M« " !) " 2 + — + U • 
The jth cumulant, j > 2, is given by 

fad - ^W)L4t* M (l+l[xJ + l) 

i=l 

fe— 1 

= ^ T yv°- l) (- + 1 ) + -^ 1) (- + -)- 

2J-i/-^ r M 2 ; 2^ v 4 2 ; 

i=0 

Differentiating ( gJBD implies ^C*- 1 )^) = ^^{z) + ^^(z + |) and therefore 

^(n, 1) = + £ ^ (1/2 + (2-20) 

8=1 

hold. The duplicity formula for 2; = | implies (f) = V* (§) — 20 (|)> where V* (z) = 

16 J2m=0 (4m+l) 2 = ^ ^m=0 ( (2m+l) 2 + (2m+l) 2 ) = 2 Z] m =0 (m+±) 2 + ^ Sm=0 (2L+I) 2 = 2V^ (l) + 

8K with Catalan's constant K, resulting in = T ~ 2ir - With and ( 1234) we 

can conclude 

r 2 (n, i) = (|) + £ (i/ 2 + = t " 2ir + log(A;) + 2 + log(2) + 1 + ' (2 ' 21) 
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For every j > 3, the first summand can be bounded using (I2.17P 



< (j - i)'(f)^ 1 + u - mq y = u - i)^(f) , 

and the remaining sum in (I2.20p is the same as in the GUE case: With (I2.13P we have 

k 

£ 1>W (1/2 + i) + (I) = _2(j - l)^- 2 \\) + (2k + ltyCi-U (1/2 + * + 1) + O(i) 



i=l 



Applying (I2.17P we obtain 



< const (j — 1)!, which implies 



the bound for the jth cumulant, j > 3. 

In the case of n = 2k even, we have to study the asymptotic behaviour of the 
hypergeometric function (see (12.51) ): F^^jl, — |; n+ ^ +s ; := 1 + Y^.=i x mi denoting 
/l+s\(m) (s\{m) 

— , 2 , r by x m . Each x m is of order 0(n~ m ) and, for s G [0,2) and n large 

/ n+i+i ^M 2 m m! v ; L ; 

enough, the hypergeometric function takes values in the interval (—1,1). Therefore we can 

study the power series of the logarithm and get 



logF 



s + 1 s n + 1 + s I 
2 '~2' 2 '2 



= logM+^a; m J = $^ m + ^(-l)'y(^0' 

^ m=l ' m=l 1=2 m=l 



We differentiate each x m via the quotient rule and the product rule in the enumerator. Setting 
s — 0, the only remaining term in the enumerator is the one where we differentiate the factor 
— |. Thus the square of the denominator cancels out. The derivative of x m equals a constant 
. It follows that the sum over I is of order 0(n~ 1 ), too. Similarly we obtain 



times 



i i 



that for every j > 1 



dP ( s + 1 s n + l + s 1 
— logFI - ;- 



0(l/n). 



Thus with O and ( 1230]) it follows that 
r x (n,l) : 



n + l_ d. + l sn + l + sl 

bg(2) + — logF , — ; ;- 

2 &v; c/s & \2'2' 2 '2 



s=0 



1 , ,lv 1 , ,71+ lv 

+ 5*( 5 )-5*(— ) 



771=1 

n + 1 



/ \ ^/1\ 1 , /1\ 1 , /Ti+1\ 1 , /1\ n n + 1, 
log(2) + O(-) + -^(-) - -'«>(—) - -«)(-) - - + — log(-) 



n 



^log(n)-| + ~log(2) + 0(l/n) 



MODERATE DEVIATIONS FOR THE DETERMINANT 



11 



and by fl2l3|) 

r,-(n,l) = — logF , — ; ;- + — V> u ^ 1} ( ) 

n ' ' dsJ V 2 2 2 2 / ,=o 2J 2 2J 2 

+ ^^'-^(1/2+771) 

m=l 

= (\) ~ j 3 ^- l) C^ 1 ) + £ ^ (l/ 2 + m) + O (1/n) • 

m=l 

Note that the only difference to the case n — 2k + 1, see (I2.20p . is the summand 1 (^r)) 

which is of order 0(l/n). Therefore the second and higher cumulant satisfy the stated bounds 
for all n. □ 

Good estimates on cumulants imply asymptotic results for the log-determinant of GUE and 
GOE ensembles, respectively Before we state our results, we remind the reader on Cramer-type 
moderate deviations and a moderate deviation principle. The classical result due to Cramer is 
the following. For independent and identically distributed random variables Xi, . . . ,X n with 
E(Xi) = and E(Jf^) = 1 such that Ee*°' Xl ' < c < oo for some to > 0, the following expansion 
for tail probabilities can be proved: 

P(W n > x) 



1 - $(x) 



l + 0(l)(l+a^)/>/n 



for < x < n l l % with W n := (X 1 + • ■ ■ + X n ) / y/n, $ the standard normal distribution function, 
and 0(1) depends on c and to- This result is sometimes called a large deviations relation. Let 
us recall the definition of a large deviation principle (LDP) due to Varadhan, see for example 
Dembo and Zeitouni [5]. A sequence of probability measures {(// n ),n G N} on a topological 
space X equipped with a a-field B is said to satisfy the LDP with speed s n /* oo and good rate 
function /(•) if the level sets {x : I(x) < a} are compact for all a G [0, oo) and for all T G B 
the lower bound 



and the upper bound 



liminf — log/i n (r) > — inf I(x) 

n->oo s n a-eint(r) 



limsup — log/i n (r) < — inf I(x) 

ci(r) 



hold. Here int(r) and cl(r) denote the interior and closure of V respectively. We say a sequence 
of random variables satisfies the LDP when the sequence of measures induced by these variables 
satisfies the LDP. Formally a moderate deviation principle is nothing else but the LDP. However, 
we will speak about a moderate deviation principle (MDP) for a sequence of random variables, 
whenever the scaling of the corresponding random variables is between that of an ordinary Law 
of Large Numbers and that of a Central Limit Theorem. 
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We consider 



as well as 



log | det*£|- 1^(71,0) 



op 



log | det Xj j\ - flog n+ f 
-Uogn 



for /3 = 1, 2 



for = 1,2. 



(2.22) 
(2.23) 



Theorem 2.3. For (3 = 1,2 we can prove: (1) Cramer-type moderate deviations: There 
exists two constants C\ and C2 depending on (3, such that the following inequalities hold true: 



log 



P(W ntP > x) 



and 



log 



1 - $(x) 



< a 



l + x s 



$(— X) 



< Co 



for all < x < C\0~$. On all cases erg is of order y/\ogn. 

(2) Berry-Esseen bounds: We obtain the following bounds: 

sup \P{W n> p <x)- $(ar)| < C(/3)(logn)- 1/2 , sup \P(W n> p < x 



$(x)| < C(/3)(\ogn)- 1/2 . 

(3) Moderate deviations principle: For any sequence (a n ) n of real numbers such that 1 <C 
a n <C o~p the sequences (~~W / n,/3) n and {^~^n,p) n satisfy a MDP with speed a\ and rate function 
I(x) = y> respectively. 

Remark 2.4. The Berry-Esseen bound implies the Central Limit Theorem stated in (II Ah . 
The statement of the central limit theorem in [H] was given differently. In section III, they 
considered a variance of order 2a 2 = 4^, meaning that the spectrum of the GUE model is 
concentrated on a finite interval (the support of the semicircular law). Then the D is the 
determinant of the rescaled (!) GUE model, given a f logn + nlog2 summand in addition to 
the expectation — n(| + log2) +0(^) they stated in [T4J (43)]. This is actually the expectation 
in (11.41) . Choosing the variance a 2 = j- in the case (3 = 2 implies that we have to rescale each 
matrix-entry Qj by Qj/ (2-y/n) and hence the determinant of the rescaled matrix is 2 n n n / 2 times 
the determinant of the matrix X 2 . 

Proof. With the bound on the cumulants (12. 7p we obtain that |rj(W nj2 )| < 7£j. With a\ > 
±(7 + 2 log 2 + 1) we get 



1^(^)1 <j\ 



1 



7-2 



a\~ 2 (7 + 2 log 2 + 1) 



<j\ 



J-2 



5< j\ 



a. 



5 \J-2 



(?2 



< 



with A = 02/5 for all n > 2. With Lemma 2.3 in [?] one obtains 



P(W n , 2 >x) 

1 - $(x) 



exp(L(x)) ( 1 + q\ip{x 



x + 1 

'"ST 
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P{W n , 2 <-x) f x + 1 



and 

for < x < Ax, where Ax = v^A/36, 

/ 60(1 + 10Afexp(-(l -x/A x ) VAT)) 

= ' 

gi,g2 are constants in the interval [—1,1] and L is a function defined in [?, Lemma 2.3, eq. 
(2.8)] satisfying |-L(x)| < for all x with \x\ < Ai. The Cramer-type moderate deviations 
follow applying (8j Lemma 6.2]. The Berry-Esseen bound follows from [?, Lemma 2.1] which is 

18 1 

sup |P(W„ )2 < x) — $(x)| < — = const 



A 1 (logn)V2- 



The same Berry-Esseen bound follows using the asymptotic behavior of the first two moments. 
Finally the MDP follows from [61 Theorem 1.1] which is a MDP for (—W nt 2) for any sequence 
(a n ) n of real numbers growing to infinity slow enough such that a n /A — > as n — > oo. Moreover 
(a _ ^ / "' 2 )n anc ^ {a~Wn,2) n are exponentially equivalent in the sense of [SJ Definition 4.2.10]: 
with W n> 2 := log l detX "|~2 lQ g"+2 we jj^g \W Ht 2 - W nt2 \ ->• as n ->• oo, and it follows that 

(a~^ /ra ' 2 )n anc ^ (o _ ^ /7n ' 2 )n are ex P onen t ially equivalent. By Taylor we have | ^-(W nt 2 — ^,2) | = 
o(l) Wn,2 and hence the result follows with |5] Theorem 4.2.13]. □ 

Next we will consider the following class of random matrices. Consider two independent 
families of i.i.d. random variables (Zij)i<i<j (complex-valued) and (li)i<i (real- valued) , zero 
mean, such that EZ^ 2 = 0,E|Z 12 | 2 = 1 and EYj 2 = 1. Consider the (Hermitian) n x n 
matrix M n with entries M*(j,i) = M n (i,j) = Z it j for i < j and M*(i,i) = M n (i,i) = Y%. 
Such a matrix is called Hermitian Wigner matrix. The GUE matrices are the special case with 
complex Gaussian random variables N(0, l)c in the upper triangular and real Gaussian random 
variables N(0, 1)r on the diagonal. 

We say that a Wigner Hermitian matrix obeys Condition (C\) for some constant C if one 
has 

E|Z M -| C < d and E\Y t \ c < C x (2.24) 

for some constant C\ independent on n. Two Wigner Hermitian matrices M n = (Ci,j)i<i,j<n 
and M' n = (CI j)i<ij<n match to order m off the diagonal and to order k on the diagonal if one 
has 

E((Re(C M )) a (Im(C M )) 6 ) = E((Re(Ci i )) a (Im(Cy) 6 ) 

for all 1 < i < j < n and natural numbers a, b > with a + b < m for i < j and a + b < k for 
i = j. 

Applying [261 Theorem 5], the Four Moment Theorem for the determinant, we are able to 
prove an MDP for the log-determinant even for a class of Wigner Hermitian matrices. For any 
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Wigner Hermitian matrix M n consider 

log | det M n | — | log n! + | log n 



W n :~- 



\ logn 



Theorem 2.5 (Universal moderate deviations principle). Let M n be a Wigner Hermitian ma- 
trix whose atom distributions are independent of n, have real and imaginary parts that are 
independent and match GUE to fourth order and obey Condition (Ci), f)2.24p . for some suffi- 
ciently large C , then for any sequence (a n ) n of real numbers such that 1< a„ < y / \ogn the 
sequence (^-W / n) n satisfies a MDP with speed a\ and rate function I(x) = If M n matches 

GOE instead of GUE, then one instead has that (^^^W n ) n satisfies the MDP with same 
speed and rate function. 

Proof. Let M n be the Wigner Hermitian matrix whose entries satisfy the conditions of the 
theorem and M' n denotes the GUE matrix. Then [26, Theorem 5] says that there exists a small 
c > such that for all G : R R + with | £jG(x)\ = O(n co ) for j = 0, . . . , 5, we have 

|E(G(log|det(M n )|)) -E(G(log|det(JlO|))| < n~ c ° 

We consider for any b, c G K the interval I n := [b n , c n ] with 



b n := ba n \j)^\ogn + ^logn! - ^logn and c n := ca n \j ^logn+ ^logn! - ^logn 

With 1+ := [b n - n~ c °/ 10 ,c n + n" Co/1 °] and I~ := [b n + n~ c °/ 10 ,c n - n -00 / 10 ] we construct a 
bump function G n : K — > M + which is equal to one on the smaller interval I~ and vanishes 
outside the larger interval J+. It follows that P(log | det(M„)| G J n ) < EG n (log | det(M n )|) 
and EG n (log | det(M^)|) < P(log | det(M^)| G J+). One can choose G n to satisfy the condition 
| -£jG n (x) | = 0(n c °) for j = 0, . . . , 5 and hence 

P(log | det(M n )| G I n ) < P(log | det(M;)| G J+) + n^ c °. (2.25) 

By the same argument we get 

P(log | det(M;)| G J") - n-<* < P(log | det(M n )| G J n ). (2.26) 

With P(^W n G [6,c]) = P(log|det(M n )| G /„). With fl2^5|) and Lemma 1.2.15] we see 
that 

lim sup — log P(W n /a n G [b, c]) < max|limsup — logP(log | det(M^)| G 7^"); limsup — logra - ' 



n— >oo 



For the first object we have 



1 1 / 1 ~ 

limsup — log P (log | det (M' n ) | G /+) = limsup — log PI — W„ )2 G [b - rj(n), c + 77(71)] 
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with r](n) := n c °/ 10 (a„y^ | logn) 1 — > as n — > oo. Since c > and log n/a 2 — )■ oo for 
n — 7- oo by assumption, applying Theorem 12.31 we have 



1 x 2 
limsup — \ogP(W n /a n G [b,c]) < - inf — . 

x&[b,c] Z 



n— too 



Applying (I2.26j) we obtain in the same manner that 



1 x 2 
limsup -5- log P(W n /a n G [b,c]) > - inf — . 
at, x&[b,c] 2 



71— >00 



The conclusion follows applying Theorem 4.1.11 and Lemma 1.2.18]. □ 

Remark 2.6. The bump function G n in the proof of Theorem 12.51 can be chosen to fulfill 
|^-G ra (x)| = 0(n c °) for j = 0, . . . , 5 uniformly in the endpoints of the interval [b, c). Hence 
the Berry-Esseen bound in Theorem 12.31 can be obtained for Wigner matrices considered in 
Theorem E2D 

sup \P(W n <x)- $(x)| < const ((logn)~ 1/2 + n~ co ). 

We omit the details. 



3. Non-symmetric and non-Hermitian Gaussian random matrices 

As already mentioned, recently Nguyen and Vu proved in [19J, that for A n be an n x n 
matrix whose entries are independent real random variables with mean zero and variance one, 
the Berry-Esseen bound 

sup \P(W n <x)- $(x)\ < \og- 1/3+o{1) n 

xeR 

with 

log(|detA n |) - ilogfn- 1)! . . 

W n := ; 2 6V (3.1) 

yi lo s n 

holds true. We will prove good bounds for the cumulants of W n in the case where the entries 
are Gaussian random variables. Therefore we will be able to prove Cramer-type moderate 
deviations and an MDP as well as a Berry-Esseen bound of order (logn) -1 / 2 (and it seems 
that one cannot have a rate of convergence better than this). In the Gaussian case, again 
the calculation of the Mellin transform is the main tool. Fortunately, the transform can be 
calculated much easier. 

Let A n be an n x n matrix whose entries are independent real or complex Gaussian random 
variables with mean zero and variance one. Denote by A\ the transpose or Hermitian conjugate 
of A n according as A n is real or complex. Then AnA^ is positive semi-definite and its eigenvalues 
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are real and non-negative. The positive square roots of the eigenvalues of A n J^ n are known as 
the singular values of A n . One has that 



J] A, 2 = det(A n A{) = | det A n \ 2 = J] 

i=l i=l 



X 



2 

i\ i 



where Aj are the singular values and Xi are the eigenvalues of A n . Now AnA^ is called Wishart 
matrix. For the real case we consider independent N(0, 1)k distributed entries, for the complex 
case we assume that the real and imaginary parts are independent and N(0, 1)r distributed 
entries. These ensembles are called Ginibre ensembles. One obtains for the joint probability 
distribution of the eigenvalues of A n A\ t on the density 



7 v 2 

ZJn ^ i=l i=l i<j 

with = 1 for the real and = 2 for the complex case and Z n ^ being the normalizing constant 
(see for example [21 Chapter 7]). As a result the Gaussian joint probability density for the 
singular values Aj gets transformed to 



1 n n n 

Q n ^K . . . , \ n ) : = —— exp (-| X i) II ^ II l A i - \ 



n ' PK ' i=l i=l i<j 



2 1/3 



with 

:=/•••/ ex P (-f E^n^^ni^-^n^ ^ 

^ 2=1 i=l i<j i=l 

Now the Mellin transform of the probability density of the determinant of A n is given by 

MnAs) = r ■ ■ ■ r iai • • • a^-u^ai, . . . , a„) n ^ = g^^+^n . 

But using the Selberg identity of the Laguerre form, [T7"t formula 17.6.5], we obtain for the 
moment generating function M n ^(s) = JA n ,p( s ~ 1) : 

JK^W = ( ? ) II r((^)/2) ■ (3 ' 3) 

This formula makes even sense for = 4, where A n is a quaternion matrix and A^ n denotes 
the dual of A n (see [T?] Section 15.4] for a discussion of the definition of a determinant in this 
case). We will concentrate on the real case (3 = 1. The results of the following theorem can be 
stated and proved similarly in the two other cases (3 = 2,4. We omit the details. We consider 
W n as in (13. ip and 



^ ;= WB|U . ^l-M^ (3 _ 4) 



log 


det A n 


- E(log 


det A n |) 


V(log 


I det A n |) 


,1/2 
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Theorem 3.1. Let A n be an n x n matrix whose entries are independent real iV(0, 1)r random 
variables. Then we have: (1) Cramer-type moderate deviations: There exists two constants 
C\ and C*2 depending on (3, such that the following inequalities hold true: 



log 



P(W n > x) 



and 



log 



1 - $0) 
P(W n < -x 



<a 



1 + x 3 



&(-x) 



<C, 



1 + X s 



for allO<x< CiV(log | detA n \)^ 2 . 

(2) Berry-Esseen bounds: We obtain the following bounds: 



sup \P(W n <x)- $(x)| < C(/3)(logn) 



-1/2 



sup \P(W n <x)- $(x)| < C(/3)(logn) 



-1/2 



(3) Moderate deviations principle: For any sequence (a n ) n of real numbers such that 1 <C 
a„ <C ap the sequences (^ _ W / n) n and (^Wn) n satisfies a MDP with speed a 2 n and rate function 



I(x) 



-, respectively. 



Proof. With (I3.3P we are able to estimate the cumulants Tj(n) of log | det A n \. The calculations 
will benefit from a few results presented in the proofs of Lemma [2TT1 and Lemma \2. 21 Therefore 
we restrict ourselves to the major steps of the proof. We denote by ip the digamma function 
and by fcGN, the polygamma function (see Lemma [2.ip . With (13. 3p we have 

n n 

r 1 (n) = ^logn + -^^(z/2) and Tj(n) = —^2ip u ~ 1) {i/2) for j > 2. 

i=l i=l 

For n = 2k + 1 we have \ YJU V>(V 2 ) = \ (E?=o ^(V 2 + i) + E?=i • Using <^M) the first 
summand is equal to -§ + f logA;+ |log & + \ijj(l/2) +0(l/k). With ^(1 + x) = ip(x) + \ (see 
Lemma I2TT1) one obtains that ip(i) = ^>(1) + E}=i j. Thus applying (12.91) we have ~ Ei = i if)(i) = 
| \og(k — 1) — | + const +0(l/k). Summarizing we get 



ri(2fc + i) 



—A; + k log /c + - log k + const +0(1/ k) 



n n 
•-(l + log2) + -log(n 



1) - - log(n - 1) + const +0(l/n). 



Therefore the leading term of the expectation of log | det A n \ is log ((n — 1)!) . In the case n 
one obtains the same order. For Tj(2k + 1) with j > 2 we proceed as following: 



2k 



1 2fc+l , k k 

r^k +i)=-x; i^cva) = - ^- i} (i/2) + E^ (i-i) (V2 + o + x 

i=l ^ i=l i=l 
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Take the representation (127121) to see that = - 1)! YlZ=i such that 

k , k 1 oo 1 \ 

E^'W = (-W-i)'(E^r + * £ sf) 

i=l \n=l m=fc+l ^ 



= -(j - 1)^ (J - 2) (1) + 0(l/ife) + k^~ l \k + 1). 
With the help of (12~T3"|) we obtain for j > 3 that 

r,H = ^l T ^- 1 )(i/2) - l(j-i)(^^ 2 )(i/2) + ^'- 2) (i)) 

+ ^(2k + 1)^^(1/2 + fc + 1) + ^k^~ l \k + 1) + 0(l/k). 

With (I2.17P we are able to bound the cumulants in a similar way as in the proof of Lemma 12.11 
and obtain | T^- (n) | < const j!. Moreover with f)2.14p we obtain for the variance 

1 1 7T 2 

r 2 (n) = - log n + - (7 + 1 + y ) + 0(l/n). 

Therefore the leading term of the variance of log | det A n \ is |logn. Now the theorem follows 
exactly as in the proof of Theorem 12.31 □ 

Remark 3.2. Let A n be an nxn matrix whose entries are independent complex and quaternion, 
respectively. Then W n and W n as defined before satisfy Cramer-type moderate deviations, 
Berry-Esseen bounds and a moderate deviations principle. This can easily be checked noting 
that, for (3 = 1,2, 4, „ 

rf(n) = ^iog(|)^ } + iE^- 1) (T 

is of order -L log(n): For j3 = 2 we have already bounded these summands in the proof above. 
In the case /3 = 4 use (12.181) and its derivatives to see, that the cumulant can be represented 
via sums of ^-^(i) and + 1/2). 

Remark 3.3 (Trace-fixed ensembles). In [15], the authors considered fixed-trace Gaussian 
random matrix ensembles (real-symmetric and Hermitian ones). Here the trace of the matrix 
is kept constant with no other restriction on the matrix elements. These ensembles are shown 
to be equivalent as far as finite moments of the matrix elements are concerned. Especially, the 
Mellin transform of the fixed-trace Gaussian matrices can be deduced from the Mellin transform 
of the Gaussian orthogonal and unitary ensemble, respectively, see [T5J formulas (17), (20) and 
(22)]. Hence it is expected that the distribution of the log-determinant of these ensembles is 
asymptotically Gaussian with a variance of order logn. We would be able to deduce the results 
in Theorem 13.11 for the Gaussian trace-fixed ensembles by the same technique. We omit the 
details. Remark, that universal limits for the eigenvalue correlation functions in the bulk of the 
spectrum for fixed trace matrix ensembles are considered in [12] , [13] . In this case, the class of 
matrices are of nondeterminantal structure. 
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